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Abstract 



Various differential cross sections of high-energy photon splitting in the electric 
fields of heavy atoms are calculated exactly in the parameter Za. The consideration 
is based on the quasiclassical approach applicable for small angles between all photon 
momenta. The expressions obtained are valid for arbitrary transverse momenta of 
final photons . The detailed investigation of the process is performed taking into 
account the effect of screening . The exact cross section turns out to be noticeably 
smaller than the result obtained in the Born approximation. 



1 Introduction 



The first successful observation of high-energy photon splitting in the electric fields of 
atoms has been recently performed in the Budker Institute of Nuclear Physics. A crystal 
of Bi 4 Ge30i2 has been used as a target. At the present time, the data processing is 
almost completed and preliminary results are published in |]]. Theoretical and experimental 
investigation of this nonlinear QED process is important as a new test of QED in strong 
external fields. It also gives a possibility to understand the role and the structure of higher 
orders of the perturbation theory with respect to the external field since, as is shown in the 
present paper, the exact in Za cross section {Z\e\ is the nucleus charge, a = e 2 /4tt = 1/137 
is the fine-structure constant, e is the electron charge) differs essentially from that obtained 
in the Born approximation, i.e. in the lowest order in Za. 

The cross section of photon splitting was found in 0, |3| in the Born approximation. In 
the same approximation, an essentially simpler form of the cross section was obtained in [|J 
with the help of the Weizsacker- Williams method providing the logarithmic accuracy. Using 
the analytical results of @, §], the cross section of the photon splitting was investigated 
numerically in || 

The Coulomb corrections represent the difference between the exact ( with exact account 
for an external field ) cross section of the process and the result obtained in the Born 
approximation. Basing on our experience concerning Delbriick scattering, we expected a 
measureable effect in the photon splitting too. In recent papers [0, ^| general formulae 
for the high-energy photon splitting amplitudes have been derived exactly in Za for small 
angles ji and fa between the momenta k2 , k3 of the final photons and the momentum 
ki of the initial one ( U{ — |kj| > m , m is the electron mass). It is the region of 
small angles, which makes the main contribution to the total cross section of the process. 
Additionally, small angles and high photon energies allow one to apply the quasiclassical 



approach, developed in [|, [UJ at the consideration of Delbriick scattering (see recent review 
PH ) . This approach is based on the use of the quasiclassical Green function for the Dirac 
equation in the electric field of an atom, which is equivalent to a summation for all orders 
of the perturbation theory with respect to the external field. The quasiclassical approach 
greatly simplifies the calculation of amplitudes. 

First theoretical results concerning Coulomb corrections in the process of photon split- 
ting were obtained in|| for large transverse momenta of the final photons as compared to 
the electron mass: |k 2 j_| = 1^2/2 3> m and |k 3 j_| = u 3 f 3 3> m. It turned out that in this 
kinematical region the Coulomb corrections lead to the significant decrease of the cross 
section. They become noticeable starting from relatively small Z and reach several tens 
per cent for heavy atoms. 

In the present paper we obtain the exact cross section of the process in a simple form, 
valid for arbitrary transverse momenta of the final photons . Using this result, we examine 
numerically the role of the Coulomb corrections in various differential cross sections. We 
discuss the case of a pure Coulomb field as well as the effect of screening. 
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2 Amplitudes of the process 



As was shown in ]7|, |j, it is convenient to present the initial expression for the photon 
splitting amplitude in the form containing the Green functions D(x,x') of the "squared" 
Dirac equation: 

D(x, x') = (x\ l/(V 2 -m 2 + iO) \x ') , 

where V = r y' J '(id fl — g^U^)), U(r) is the potential energy of an electron in an external 
field, 7^ are the Dirac matrices. Then the amplitude M is splitted into a sum of the 
perturbation-theory diagrams, containing either three or two Green functions D(x, x')\ 
M = M (3) + M^. The term Af® is given by 

M (3) = -e 3 J — J dridr 2 dr 3 exp[i(kiri - k 2 r 2 - k 3 r 3 )] x (1) 
Trj [{-eA - 2eip)£)(n, r 2 \e - cu 2 )] [{efk 2 - 2e 2 p) J D(r 2 , r 3 |e)] x 
[(e* 3 k 3 - 2e 3 p)D(r 3 , n \e + u 3 )] \ + (k% <-> k% , e 2 <-> e 3 ) . 



Here ei and e 2 3 are the polarization vectors of the initial and final photons, e = = 
— erf, and the operator p = — «V differentiates the Green function D with respect to its 
first argument. The term reads 



j^J dY 1 dr 2 Tr|exp[i(kiri - k 2 r 2 - k 3 r 2 )] e^e 3 x 



(2) 



[(-eiAn - 2eip)D(r ls r 2 \e - wi)] D(r 2 , r x |e)] + 



exp[z(kiri - k 2 r 2 - k 3 ri)] eie 3 x 



D(r u r 2 \e - u 2 )\ [{e*js 2 - 2e* 2 p)D{r 2 , v 1 \e)\ + ^ k% , e 2 e 3 ) 



As was pointed out in 0, the effect of screening is important only for the lowest (Born) 
approximation in Za. Therefore, we start from the case of a pure Coulomb potential 
U(r) = —Za/r and then multiply the Born contribution to the amplitude by the atomic 
form factor to take the effect of screening into account. 

It is convenient to perform the calculations in terms of the helicity amplitudes 
MxiA 2 A 3 (ki, k 2 , k 3 ). We direct the z axis along ki and introduce the vectors f 2 = k 2 ±/u 2 
and f 3 = k 3 _|_/co> 3 (|f 2)3 | 1). The z component of the polarization vectors ej can be elim- 
inated owing to the relation e^kj = which leads to e z = — e_i_k_i_/cG>. After that within the 
small-angle approximation one can neglect the difference between the vectors (e 2 ^)± and 
the polarization vectors of photons, propagating along the z axis and having the same helic- 
ities. Therefore, the amplitudes M\ 1 \ 2 \ 3 (ki, k 2 , k 3 ) contain only the transverse polarization 
vectors e and e*, corresponding to the positive and negative helicities, respectively. It is 

sufficient to calculate three amplitudes, for instance, M_| (ki,k 2 ,k 3 ) , M+++(ki, k 2 , k 3 ) 

and M++_(ki, k 2 , k 3 ). Other amplitudes can be obtained by the substitution 

M + _ + (kx, k 2 , kg) = M ++ _(kx, k 3 , k 2 ) , M_ A2A3 (k!, k 2 , k 3 ) = M +A2 A 3 (k 1 , k 2 , k 3 ) (e <-> e*) , 

where A denotes the helicity opposite to A. 
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It was shown in that the main contribution to the term is given by the region 
where Z\ < and at least one of z 2 and z 3 is positive. Similarly, the main contribution to 
the term is given by the region z\ < and z 2 > 0. General formulae for these terms 
at arbitrary relations between the electron mass and the transverse momenta of the final 
photons are presented by Eqs. (9-14) and Eqs. (16-17) in 0. Performing an integration by 
parts in the expression for the term M^ 3 ' as it has been done when deriving Eqs. (18-19) in 
0, but keeping now the terms containing the electron mass, we obtain for the amplitude 
M of the photon splitting 

M = Mi + M 2 + 5M , (3) 

where Mi corresponds to the contribution to from the region Z\ < < z 2 , z 3 , while 

M 2 corresponds to that from two regions Z\ < z 3 < < z 2 and z\ < z 2 < < z 3 . The 

quantity SM is a sum of the term and the integrated terms arising when is 

integrated by parts . 



For the term Mi, we have 

. , Ll>2 oo oo L 

lids P C C f Cli -t\j 

(Mx) Xl x 2 x 3 = — -= / ek 2 k z de / gLRi / dR 2 / — 4 x (4) 

o o o o 

/ \ 2iZa 

J J dq 2 dq 3 T AlA2A3 e**Im ( + (u 2 uj 3 , k 2 k 3 , A 2 A 3 ), 

where q 23 are two-dimensional vectors lying in the xy plane, k 2 = cj 2 — e , K3 = uj 3 + e, 
L = R 2 uj 3 k 2 /uj 2 k 3 , 







[(5- 





2 2i? Wl U 



u^R y ^ ' 2 

R = R 2 -R 3 , f 23 = f 2 -f 3 , Q = q 2 + q 3 , A = ^ 2 f 2 + ^ 3 f3- 

The function T for different polarizations is: 

T+-- = ^(eQ)(eQ 2 )(eQ 3 ), (6) 



— + — ( e Q)[( e Q 2 )( e Q 3 ) + e ~^Qa + -ftQ 



R\R? \ e k 3 J ~ k 2 \k 3 R eR\ 

where Q 2 = Q + £-R 2 f 2 3 and Q 3 = Q + eR 3 i 23 . 
The term M 2 is given by 
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oo oo 



(M 2 ) Al 



J de J dRi J dR,2 J 



dR 



dq 2 dq 3 T Xl x 2 x 3 e i$ Im f — 





2iZa 







RoR 



ix 



(7) 



+ (u>2 <-> co> 3 , k 2 <->• k 3 , A 2 A 3 



where L = RxUJ^^j oj\E , R = R\ + R3 , and the functions $ and T can be obtained from 
$ and T in ©,(0) by the substitutions 



^2,3 -* -Q2,3 > w l ^ ^2 , ^3 -> ~^3 , ^3 «"> £ , 
R\ ^> R2 , -R3 — > — -R3 > f23 ^ — f3 j f 2 — ► — f 2 , 



so that 

rri rri rji rji rji 

-H ~~ * -H > ■*■+-+ ~* ■*■+-+ 7 J ++ 







22 



2iZa 



where 





+ = (eQ) 


I 






—U>3 








= (e*Q) 


/ de—^- 









<"' ! - /tfe-^e^ 2 



L 3j , 



F++- = F + - + (UJ 2 <-> ^3 , f 2 ^ f3) 



1 1 Q ^2^3^2^3.2 D 

1 1 Q 2 WiW 3 K 2 £ f 2 p 



(8) 



> T++ _( e ^e*), T ++ _^T +++ (e^e*). 
For the last contribution to the amplitude M (|3|) we obtain 

(5M) Xl x 2 \ 3 = J dRi J J J rfq 2 dq 3 F AlAaA3 Im 



(9) 



(10) 



(K 2 q 2 - K 3 q 3 , A) m 2 

— (-Ki + #2) 



■/t 2 q 2 - gq 3; A) 

lJ 2 



777 



The expressions derived are still rather complicated and require further transformations. 
Let us pass from the variables q 2 , q 3 to Q = q 2 + q 3 and q = q 2 — q 3 and use the identity 
(see Eqs. (22,23) in 0) 




2iZa , . /i * 1 \ 2iZa 

./|> exp( _i qQ)to g±^ . 



Additionally, the parametrization 



4 



6Xp( ^2~k' ) =iRl I^ exp ( _i ~f~ ~~ i( ^ x ) ' ( 12 ) 
where x is a two-dimensional vector, is used to calculate the term M\. 

After that it is easy to take first the integrals over Ri, Q, and then over i? 2 and _R 3 . 
The calculation of the terms M 2 and 6M is carried out analogously. Performing, finally, 
the shift x — > x — q/2, we obtain only two kinds of integrals over q to be taken: 



G = m 



cos r 



/ q+ _ q- 

x+A/2 x-A/2 



2iZa 



212 



Re 



Gi = 2m 2 



T\ + sin r 



q/2) 

'x + A/2 x-A/2 N 



(13) 



r ^q 


f e q+ 


eq_\ 


/ 2tt 




<?- ) 



e*(x-q/2) 



^2, 



-Re 



m cos r 



JF X + sin r 1 — 



2j2 
2 m 2 



9+ 

1^2 



2/Za 



where the notation 



q ± = q ± A , c± = m 2 + (x ± A/2) 2 , a = 

is introduced. The functions T\ and ^2 are given by 

dx 



\ 



cos (2Za arcsinh(aa;)) 



7o (1 

= / sin (2Za arcsinh(ax)) , = 

7o (i + x 2 ) 3/2 1 1 )J VTTW 



(14) 



We explain in Appendix I how the integrals in (|13l) have been actually taken. The represen- 
tation (p~3|) is very useful since a dependence of the functions T\ and JF 2 on x is expressed 
by the single variable a. This greatly simplifies numerical calculations. 



Finally, we obtain for the helicity amplitudes of photon splitting 



M+__ = N I dx (eG) / ^k 2 



A 



ea 



e acq 
+ 



+ 



(15) 



23/ 



k 3 (ef 



g(ef 3 ) 

(e*f 23 ) w 2 £> 3 (e*f 3 



Nldxl - A 



(eG) (k? 



(ef; 



+ 



U2 <-> cj 3 
f 2 <-> f3 

m 2 UJ 3 K2 (e*f ; 



23/ 



23/ 



+ 



^i^KsGie* (a + f 2 3Cj 2 K 3 M) u; 3 K 2 (e*G) 

H — [£CU 3 y4 + 



2?i 



^2^3 
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2e (k\ + «§) (e*f 3 ) (ea) + m 2 (2eK 2 - LOiUJ 2 )\ 



LU1UJ3K2G1 (e*c) 



V 3 



co 2 <-> 
f 2 <-> f 3 



++- 



iV / rfx 



^2 de 
2A 



k 2 u 3 (eg) 
wiX>i 



[«3 (« 2 - e) A- 



2/« 3 («| + £ 2 ) (ef 2 3) (e*a) + m 2 {u)\U) 2 - 2k 2 k 3 )] + 



k 2 uj 2 uj 3 Gi (eb) 

^1 



+ 



+ e 2 ) (e*G) ( m 2 uj 3 K 2 (ei 3 )\ u 2 u 3 K 2 Gie (a + f 3 uxe/uj 2 ) ' 

— ea H — + 



(e*f 3 



, u; 2 k 3 
de — — 
-U13 2B 



to 2 V 3 



(eg) 



[- («1 + e« 3 ) 5 + 2«a (4 + e 2 ) (e*f 23 ) (eb) + 



9 / WoGi (eb) feG) r 

m 2 (c^ 2 - 2k 2 k 3 )} + j- V ; + ^— ^[-c^fi + 

Z^l ^ 3 i^ 2 

2k 3 («| + £ 2 ) (e*f 2 ) (eb) + m 2 (cu 2 cu 3 - 2eK 3 )] + ^ ' ' ' ^ ' 



where the following notation is used 



a = x-A/2 + /t 2 f 2 , b = x + A/2 - K 3 f 3 , c = x + A/2 - rf 23 , (16) 

8e 3 Za . 22 ^ 2,2 

^ = ~~ oTo » A = m + a , 5 = m 2 + b 2 , 

V 2u; 3 / uj 3 
-n f , K 2 - e \ 2 WiW3/f 2 e j 



As was pointed out in [|7J , the quantity A 2 in the coefficient N should be interpreted as 
the squared total momentum transfer: 

A 2 = (k 2 + k 3 - k x ) 2 = A 2 + A 2 = (k 2± + k 3± ) 2 + \ ( ^t! + 2 . (18) 

4 \ u 2 uj 3 ) 

Since the functions G and G\ are independent of the energy e, the integrands in (|T5|) 
are rational functions of e, where all the denominators are quadratic forms of this vari- 
able. Therefore, the integrals over e can be expressed via elementary functions. Resulting 
formulae being rather cumbersome are not presented here explicitly ( some details of this 
integration are described in Appendix II). Performing the integration over e in fll5|) ,we 
obtain in fact a twofold integral over x for the amplitude of photon splitting since for a 
given Za, the functions T\ and T 2 (see (0)) can be tabulated separately. In the limit 
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m — > the amplitudes (|15|) coincides with those obtained previously in ||. Additionally, 
it has been checked numerically that in the limit Za — > our results ( |I5D agree with those, 
obtained in || |6| in the Born approximation. 



3 Cross section 



In the small-angle approximation ( j F 2 [ , |f 3 | <C 1) the cross section of photon splitting has 
the form 

da = |M| 2 f ^ dk3± r , (19) 
1 1 2 8 n 5 ujx(l - x) ' V ; 

where x = u 2 /uji, so that u 3 = u>i(l — x). As was mentioned above, in the general case 
of a screened Coulomb potential the lowest in Za (Born) part of the amplitude should be 
multiplied by the atomic form factor. The Moliere representation [12j for this form factor 
reads 



1-^(A 2 ) = A 2 E Z ^, (20) 

where 

ai = 0.1 , a 2 = 0.55 , a 3 = 0.35 , # = , (21) 

6i = 6 , b 2 = 1.2 , b 3 = 0.3 , fa = mZ l/3 /121 . 

To illustrate a magnitude of the Coulomb corrections, the exact and Born differential cross 
sections da / dxdk 2 ±dk 3 ± are plotted in Fig.l depending on \k 2 ±\/m for |k 3 j_| = 2m and the 
azimuth angle (the angle between the vectors k2j_ and k3j_) = 0, n. The calculations were 
performed for a screened Coulomb potential at x = 0.1, oj\ = lGeV. The value Z = 83 
(bismuth) was chosen since bismuth atoms determine the cross section of photon splitting in 
the experiment JIJ. A wide peak for azimuth angle = tt is due to small momentum transfer 
A. There is a narrow notch in the middle of this peak ( at \k2±\/m = 2 ) where the condition 
Aj_ = k2j_ + k3_i_ = is fulfilled. The width of the notch is about max(A z /m, j3o/m). Recall 
that A 2 is the longitudinal component of the momentum transfer defined by (0), and (3q 



(pi|) characterizes the effect of screening. In our example (3q is larger than A z , so the width 
of the notch is roughly flo/m = 3.6 • 10 -2 . Let us note that for u)\ ^> m the differential cross 
section, expressed in terms of k 2 j_, k 3 j_, x and ui\ depends on the energy ui\ only via A z . 
Due to this, the differential cross section is independent of uj\ outside the notch vicinity. 
The behavior of the cross section at small A is determined by the Born amplitude which is 
proportional to 1/A in this case. That is why the exact and Born cross sections coincide 
within peak region. Outside this region the Coulomb corrections essentially modify the 
cross section. The points of discontinuous slope on the curves in Fig.l are related to the 
threshold conditions for real electron-positron pair production by two photons with the 
momenta k 2 and k 3 : 

(k 2 + k 3 f = c^ 3 f 2 3 = 4m 2 . (22) 
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In Figs. 2-4 the differential cross section m 2 a Q l da/dx dk 2 ± is shown depending on 
k 2 ^/m for a screened Coulomb potential at ui/m = 1000, Z = 83 and different x, 

aa = ^l = 0.782- 10-* b. 

Solid curves represent the exact cross sections, and the dashed curves give the Born re- 
sults. The cross section exhibits a thresholdlike behavior in the vicinity of the point 
k 2 ± = kth = 2^x(l — x)m, where both conditions A = k2j_ + k3j_ = and (j22"D hold. 
Under these conditions the peak in the cross section da / dxdk 2 j_d]i3j_ seats on the bound- 
ary of the kinematic region where real electron-positron pair production by two photons 
with the momenta k 2 and k 3 is possible. The cross sections da/dxd\t 2 j_ drop rapidly for 
k 2 ± ^> m (oc l/A;!). The dotted curves in Figs. 2-4 show the difference between the Born 
and exact cross sections, i.e. they give the Coulomb corrections taken with the oppo- 
site sign. Again, as is seen from Figs. 2-4, the Coulomb corrections to the cross section 
integrated over k 3 x noticeably diminish the magnitude of the cross section. Above the 
threshold (k 2 ± > kth) this difference reaches tens per cent while below the threshold the 
exact cross section is several times smaller than the Born one. It can be explained as fol- 
lows. Above the threshold the main contribution to the cross section da/dx rfk 2 j_ is given 
by the integration region where m&x(kl ± / lui, f3 ) C A < k 2 ±. As a result, the Born cross 
section is logarithmically amplified as compared to the Coulomb corrections. Far below the 
threshold where k 2 ± <C kth, it follows from the condition A <C k 2 ± that k 3 ± rs k 2 j_ <C m, 
and the amplitude is suppressed as a power of fc| ± /m 2 . Therefore, below the threshold the 
main contribution to the cross section da/dxdh 2 ± is given by the region k%± ~ m, where 
the exact in Za amplitude drastically differs from the Born one. It is seen from Figs. 3,4 
that in accordance with previous discussion a position of the peak is the same for x = 0.2 
and x = 0.8. Nevertheless, the magnitudes of these two cross sections are significantly 
different , especially below the threshold. The explanation is the following. For the region 
/c 3 _i_ ~ m, which makes the main contribution to the cross section da/dxd\<. 2 ± below the 
threshold, and k 2 ± <C m the invariant s = (k 2 + k 3 ) 2 ~ m 2 x/(l — x). So, s <C m 2 at x <C 1 
which leads to the suppression of the cross section below the threshold (see Fig. 3). This 
is not the case for 1 -1 < 1 (see Fig. 5). However, the cross section da/dx for x = 0.2 
should coincide with that for x = 0.8 and this was checked numerically. 

Let us consider now the magnitude of the Coulomb corrections to the cross section da/dx 
integrated over the transverse momenta of both final photons. The main contribution to 
this cross section is given by the region where |k2j_|, |k3j_| ~ m. The Born contribution 
to da/dx contains large logarithm resulting from the integration over small momentum 
transfer region max(/3 , m 2 /u)i) C A « m. For j3 3> m 2 /ui the cross section da/dx is 
independent of tJi, while for /3 m 2 /u 1 it slowly grows (as lnt^i/m) when oj\ increases. 
Since the Coulomb corrections to da/dx are determined by the region of momentum transfer 
A ~ m, they do not depend on u>i for uoi ^> m . They also are insensitive to the effect 
of screening. In Fig. 5 the exact (solid curve) and the Born (dashed curve) cross sections 
a^ da/dx are plotted as functions of x for toi/m = 1000, and Z = 83. As it should be, 
the curves are symmetric with respect to the replacement x — ► 1 — x. Dotted curve shows 
the Coulomb corrections taken with the opposite sign. Note that their dependence on x is 
very weak. If x — > or x — > 1 then the cross section da/dx increases rapidly. However, the 
cross section da/dx should vanish at x = and x = 1 due to the gauge invariance of QED. 
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Actually, the cross section da/dx begins to decrease very close to the x interval end points 
(5x ~ m 2 /ouf). Therefore, the contribution from these x-range to the total cross section a 
is negligible. In our example ( Z = 83, u)\jm = 1000 ), the exact result for a is 3.9 ■ 10~ 4 b 
while the Born approximation gives 4.8 • 10 _4 b, the difference being 23%. 

In Fig. 6 the Coulomb corrections da c /dx divided by <To are shown as a function of Z 
for x = 0.7. Since their dependence on x is rather weak ( see Fig. 5), this curve allows 
one to estimate the magnitude of the Coulomb corrections for any x. It can be seen from 
Fig. 6 that starting from Z « 30 the dependence of a ~ 1 da c /dx on Z is almost linear and 
this quantity is not described by its lowest in Za approximation (oc (Za) 2 ), so that higher 
order terms are important. 

Thus, the process of photon splitting can be adequately described only with the Coulomb 
corrections taken into account. At large Z their contribution is always essential, though the 
magnitude of the Coulomb corrections depends on the type of the cross section, kinematic 
conditions, and Z. Indeed, the predictions based on the exact in Za cross section are 
in agreement with the preliminary experimental data of Jl|, while the Born results are 
noticeably different. 



r 


3± 




/ 2vr 


u 
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Appendix I 

In this Appendix we describe the calculation of the integrals in ([13]). Let us consider the 
integral 

/ \ 2iZa 

* + (x- q /2)^ (t J - < 23 > 
where q ± = q ± A. To transform this integral, we multiply the integrand in fl23|) by 

l = Jdy8{y- -^^s) = (q 2 + A 2 ) / ^ *((q - A/y) 2 ~ ±\W ~ 1)) , 

change the order of integration over q and y and make the shift q — > q + A/y. After that 
the integral over q becomes trivial, and the integral over the angle of q can be easily taken 
by means of the residue technique. As a result, we obtain 

G = -m 2 A 2 | i ^Re(^^J (yx-A/2)x (24) 
m 2 A 2 + (xA) 2 - 2xACy + (C 2 - m 2 A 2 ) y 2 ^ " 2 



C = m 2 + x 2 + A 2 /4 . 

Let us perform the substitution y = tanhs and then make the shift s —>■ s + (1/2) In c + /c_, 
where c± = m 2 + (x ± A/2) 2 . Finally, we come to the expression for G in fll3|) with T\ 
and JF 2 in the form of 



Ti = a 2 [ 
Jo 



00 , cosh s cos (2Zas 
ds 



( sinh s + a 2 



3/2 



_ • / , sinh s sin (2Zas 
j- 2 = a \ as 



2 j_ 2\ 3/2 

s + a 2 ) 



(sinh 

where a 2 = m 2 A 2 /c + c_. Making in these formulae the substitution sinhs = ax, we obtain 
the form fllij) . The quantity G± can be transformed in the same way. 



Appendix II 

A general form of the integral over the energy e reads 

de-%= , (25) 
o AB 

where 

A = A 2 e 2 + Ate + 1, B =B 2 e 2 + B x e + \ , 

and Q is the fifth-degree polinomial in e. To perform the integration, it is convenient to 
use the representation 

5 n4 ( e A 1 B 2 + A 2 B 1 \ 

AB A 2 B 2 n5 \A 2 B 2 (A 2 B 2 ) 2 J 

(n) (n) (n) (n) 

Pi g + Po _ Qi g + go 
A B 

where the coefficients pf 1 and satisfy the recurrence relations 



(n-l) 



in) (n-1) At ( n -l) (n) Pi 

Pi = Po - j-Pi > Po = 

(n) (n-l) -Dl (n-l) (n) Q\ 

(R = % ~ -b-Qx , % = e— ' 

D 2 D 2 



(o) _ A 2 V (o) _ AtV- A 2 (A 2 -B 2 



Pi = Po 



S 2 F (0) S X V - B 2 (A 2 - B 2 ) 



(0) _ ^ (0) 

9i — j w 



y = A 2 fix - AiB 2 , iy = (Ax - soy - (A 2 - fi 2 ) 2 

Thus, the calculation of the integral (^) reduces to the calculation of standard integral 

Jo d£ A 2 e* + Ate + l ' (26) 
expressed via the elementary functions. 

This is the algorithm used in our numerical calculations. 
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Figure captions 



Fig. 1. Differential cross section da / dxdk 2 ±dk 3 ± vs |k 2 j_|/m in a screened Coulomb 
potential for different azimuth angle <fi between vectors k 2 ± and k%±; Z = 83, x = 0.1, 
ui = lGeV, k 3 ± = 2m. The dashed curve (Born approximation) and the solid curve (exact 
cross section) correspond to <fi = n. The dash-dotted curve (Born approximation) and the 
dotted curve (exact cross section) correspond to = 0. 

Fig. 2. m 2 <jQ l d<7 / dxd\t2±_ vs |k 2 _i_]/m for a screened Coulomb potential, oj\/m = 1000, 
x = 0.5, Z = 83, <T is given in the text. The dashed curve corresponds to the Born 
approximation, the solid curve gives the exact result, and the dotted curve shows the 
difference between the Born cross section and the exact one. 

Fig. 3. Same as Fig. 2 but for x = 0.2 . 

Fig. 4. Same as Fig. 2 but for x = 0.8 . 

Fig. 5. The dependence of a^da/dx on x for a screened Coulomb potential, oj\/m = 
1000, Z = 83. The dashed curve corresponds to the Born approximation, the solid one 
gives the exact result, and the dotted curve shows the difference between the Born cross 
section and the exact one. 

Fig. 6. The dependence of the Coulomb corrections a^dac/dx on Z for x = 0.7. 
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